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Abstract 

Motivated by the renewed interest in studying the pion production on 
nuclei with protons at few GeV incident energies, we investigate the pion 
production in proton-proton colhsions over an energy range of 300 MeV 
to 2 GeV . Starting from a reahstic one-boson exchange model with pa- 
rameters fitted to the amplitudes of the elastic nucleon-nucleon scattering, 
we perform fully covariant calculations for the total, double and triple 
differential cross-sections of the p{p,mr'^)p and p{p,p7r^)p reactions. The 
calculations incorporate the exchange of vr, u; and a mesons and treat 
nucleon and delta isobar as intermediate states. We obtain a reasonably 
good agreement with the experimental data in the entire range of beam 
energies. The form of the covariant delta propagator, the cut-off parameter 
for the ttNN and vrA^A vertex form factors and the energy dependence of 
the delta isobar decay width is investigated. 
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1 Introduction 



The study of pion production in two nucleon collisions dates back to the mid 
fifties when one of the earliest measurements of the p{p, n'K^)p reaction was per- 
formed This was followed by several measurements of this reaction in 1960's 
and early seventies . These experiments were usually conducted in emulsions or 
bubble chambers, consequently the data invariably had poor statistics. However, 
with the advent of accelerators capable of producing intense high quality beams 
of protons with energies up to a few GeV and sophisticated detecting systems, 
it became possible to obtain good quality data on the charged as well as neutral 
pion production in proton- nucleon collisions at several beam energies [§, 0. Even 
the complete kinematical measurements of the three particles in the final state 
can now be performed [§, 0, 0- Furthermore, with improved polarisation tech- 
niques for beams and targets, the spin observables (asymmetries) are also being 
measured with great accuracy |^, 0. These measurements provide a very strin- 
gent test of the theoretical models of the pion production which is the dominant 
inelastic channel in nucleon-nucleon (NN) collisions. 

A thorough investigation of the pion production in the elementary NN colli- 
sions within a fully relativistic model is essential in the context of recent efforts to 
develop a covariant two nucleon model (TNM) to describe the A{p, tt)B reactions 

|10|, where several parameters have to be fixed from the study of the former 
reaction. The precise information about the pion production cross sections in 
elementary NN collisions is also important for the description of the dynamics of 
the heavy ion collisions within the kinetic theories (eg. BUU) ||TT], |12|. 

Although several models have been proposed to study the A^A^ NNtt reac- 
tions [|I3], only a few have taken the fully relativistic Feynman diagram approach 
0, |16|, |1^]. However, a systematic study of all the observables (total, single, 
double, and triple differential cross sections and asymmetries) of p{p,mT^)p and 
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p{p, 7i^)p reactions over a wide range of beam energies within the covariant model 
is still lacking, and the form of several ingredients of the model remains to be 
far from being determined unambiguously. Moreover, the role of a and u mesons 
have not been explicitly studied in such theories so far. 

Fully relativistic calculations could be necessary even at the beam energies 
closer to the NNtt threshold, where the non-relativistic models underpredict the 



experimental data by a factor of approximately 5 |T8|. One of the consequences 
of the relativistic effects is that they lead to an enhancement of the axial vector 
current of the NN system which can increase the pion production cross 
sections, and possibly contribute to an explanation of the above discrepancy (see 
e.g Horowitz et al. |jl3|]) to a great extent. 

The aim of this paper is to perform a detailed investigation of the p{p, mi^)p 
and p{p,piT^)p reactions using a fully relativistic Feynman diagram technique. 
We carry out our calculations within an effective one-boson exchange (OBE) 
nucleon-nucleon (N-N) scattering mechanism, which includes both nucleon and 
delta isobar excitations in the intermediate states. We consider the exchange of 
vr, o", p, and u mesons. Most of the parameters of the OBE model are determined 
by fitting to the N-N scattering data [|r^. Different delta isobar propagators 
proposed in the literature pp] , pi] , p2| have been examined in order to remove 
the confusion about their most appropriate form. The relative importance of the 
contributions of different exchanged meson in the energy range 300 MeV to 2.0 
GeV has been discussed. 

In the next section, we give the details of our model and derive the expressions 
of various amplitudes. The comparison of our calculations with the experimental 
data and the discussion of the results are presented in section 3. The summary 
and conclusions of our work are presented in section 4. 
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2 The amplitudes of pion production in One Bo- 
son exchange model 

To obtain the desired production cross section for the NN — > A'"A'"7r processes, 
we use a fully covariant method based on an effective one-boson-exchange model 
which describes at the same time the elastic nucleon-nucleon scattering. Pions 
are simply produced from the external nucleon lines (see Figs, la - Id) or if 
possible, also from the internal meson lines (Figs, le-lh). The pion production via 
formation, propagation and the subsequent decay of the delta isobar is included. 
It may be noted that in Fig. 1 we have not shown the so called 'pre-emisson' 
diagrams. However, in actual calculations their contributions are also included. 

2.1 Model Lagrangian 

The nucleon-nucleon interaction is described by the exchange of tt, a, p and 
uu mesons in term of which we parametrize the T-matrix. The corresponding 
Lagrangian densities are given by 



JC.nn = -^^Nl5l,r-id^^)^M. (1) 

CpNN = -gp^N[l^. + i^^Jt^ud^^r ■ pi'-^r^. (2) 

C^NN = -y.*iv(7M+2^<^M-^')^''*iv. (3) 

C^aNN = Qa^NT'^N- (4) 



It may be noted that we have used a pseudovector coupling for the ttNN 
vertex. In addition the couplings of the mesons to the delta resonance are needed. 
Due to the isospin conservation only p and tt mesons couple to the delta resonance. 
The Lagrangian densities for these processes are 
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C^NA = ^^,T-d^^^^r^ + h.c.. (5) 
CpNA = i f^T(9V^-9V'^)7,75^jv + h.c.. (6) 

To insure the gauge invariance of the p — N amphtude we must also include 
the diagrams le - Ih, which depict the processes where the p meson decays into 
two pions in flight. The corresponding Lagrangian density is given by 



Cp^. = X $] -p^. (7) 

In Eqs.( 1-7), \E'jv represents the Dirac spinor for the nucleon in the spin-isospin 
space. \E'^ is the coresponding Rarita-Schwinger spinor for the delta isobar. p, 
a and u represent the pion, rho, sigma and omega meson fields, respectively. T 
and T represent the isospin operator for the transition A — > A^vr and the isospin 



Pauli matrices, respectively. We use the conventions and notations of Ref. ^3 
for definitions of spinors, operators and isospin matrices, respectively. 

Since we use the Lagrangians (1-4) to directly model the T-matrix, we have 
also included a nucleon-nucleon-axial vector-isovector vertex, with the Lagrangian 
density given by 

In Eq. (8) A represents the axial vector meson field. This term is introduced 
because in the limit of large axial meson masses (m^i) it cures the unphysical 
behaviour in the angular distribution of NN scattering caused by the contact 
term in the one-pion-exchange amplitude ^Tj , if qa is chosen to be 



1 f u 

4 = 

with very large m^v) mA- 
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gA = -ruA — ■ (9) 
3 



The determination of other couphng constants used in this work in discussed 
in the next section. 

2.2 Coupling constants 

The couphng constants appearing in Eqs. (1-4) were determined by fitting 
to the elastic proton-proton and proton-neutron scattering data. In the fitting 
procedure we also take into account the finite size of the nucleons by introducing 
the form-factors 



'A? - m^' 



.i = 7T,p,a,u, (10) 



at each interaction vertex, where q is the four momentum and m the mass of the 
exchanged meson. The form factor supresses the contributions of high momenta 
and the parameter A, which governs the range of suppression, can be directly 
related to the hadron size. Since the data in the entire range of beam energies 
cannot be reproduced with an energy independent set of parameters, we have 
used the following energy dependence for the coupling constants 

giVs) = goexp{-l^/s). (11) 



The parameters {qq, A, and i) were determined [|T7| by fitting to the relevant 



proton-proton and proton-neutron data at three beam energies of 1.73 GeV ^ 2.24 
GeV , and 3.18 GeV , where a good fit to the elastic scattering amplitudes were 
obtained. Table 1 shows the values of the parameters obtained by this procedure. 
It may be noted that for the case of pion we have shown in this table the constant 

which is related to /tt of Eq. (1) as g.„ = (/^/m7r)2m7v. 

The value of the coupling constant /* for the ttA^A vertex has been determined 
from the A N + 7c decay and its value is 2.13. As in the case of NN scattering, 
we use form factors for the N — A vertices as well, which have a dipole form 



0,0 

^' = [w^) '''' 

The parameters g*, A* and A* were determined by fitting to mass differential 
cross sections for the reaction A^ + A^— i>iV + Ain the energy range of 1-2.5 GeV. 
Their values are, 



a: = 1.421 GeV, 



a; = 2.273 Gel/, 



g; = 7.4. (13) 



For the graphs le - Ih, we have taken gp^^^^ = 2gp p4| , |25|| . In this way the 
NN meson and A^A vertices are determined rather reliably. We then assume 
that their off-shell dependence is detemined solely by the form factors (10) and 
(12). 

For the delta isobar propagator we used the form given by Benmerrouche et 



al. 22 



p 
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Similar expression for the A propagator has been derived also in Ref. p6[. It 
should be noted that the form of the propagator for an interacting delta isobar 
remains the same as in Eq. (14) [^. In appendix B, we have presented a review 
of the derivation of this propagator and have discussed the difficulties associated 
with A isobar propagators presented by Williams and other authors [2^, |21|| . 

The mass of the delta isobar rriA appearing in the denominator term [p^ — m^) 
is modified by adding to it an imaginary width, which is a function of the pion- 
nucleon center of mass momentum p^. This is due to the fact that A isobar 



is not a stable particle and it decays with a width which varies with its mass 



23] • There are several prescriptions for the mass dependent delta isobar width 
given in the literature [1^, |2^, We have performed calculations by using the 
following forms 



where 



P. = • (16) 

In Eq. (16), pi is the four-momentum of the intermediate delta isobar (see Fig. 
1). used in Eq. (15) is obtained from Eq. (16) by substituting pf = m\. 
The constant Fq is the free delta width; its value is taken to be 0.120 GeV and 
e = 0.16. This form has been used by Dmitriev et al. Verwest |TB[ has 

used a somewhat different form which is given by 



p'^ \ I ^p'^ + ml + mN 



r«.(p;) = r„ ( i|j 11 V--- 1 , (17) 



where Fq has the same value as given above. 

In the next subsection we present the expressions for the amplitudes of various 
graphs as shown in Fig. 1. 

2.3 Amplitudes and cross sections 

After having established the interaction Lagrangians and the form of the delta 
isobar propagator, we can now proceed to calculate the amplitudes corresponding 
to the various diagrams associated with the reactions p{p,mT^)p and p{p,piT^)p. 
The Feynman rules for writing down these amplitudes are well known (see e.g. 
[Q). The isospin part is treated separately which gives rise to a constant factor 
for each graph. In table 2 we give the values of these factors for the post-emission 
(shown in Fig. 1) as well as for the pre-emission graphs corresponding to the 
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emission of 7r+ and 7r° mesons. Values for the pion emission from the intermediate 
states (Figs, le-lh) are shown in table 3. It may be noted that in these tables 
the values for the tt"*" case are calculated by assuming particles 3 and 4 as neutron 
and proton respectively. Isovector corresponds to the exchange of tt and p mesons 
while isoscalar to that of a and ui mesons. 

In the following we give the expressions for the amplitudes for the pion pro- 
duction processes via excitation of both nucleon and delta isobar intermediate 
states for one graph eg. Fig. Ic. We also give expression for one diagramm (Fig. 
If) where pions are emitted from an intermediate meson. The amplitudes for 
other similar graphs can be written in a straight forward manner. 

(i) Nucleon intermediate state and pion exchange 



An{c) = -QUc) (—] ^(PsHlxqMPiWq) 

X'0(p4)757X^iv(Pi)757r^?XP2), (18) 

where DT^{q) and DN{pi) are the propagators for the exchanged pion and the 
intermediate nucleon respectively, which are defined as 

DM) = -r^^ (19) 

n ( \ ■ Pirr/'' + .„„s 

DNipi) = I , • (20) 

Pi '''■N 

Various momenta appearing in Eq. (18) are defined in Fig. 1. The intermediate 
momenta are given hy q — pi - and Pi — Ptv + Pi- V' is the Dirac spinor in the 
spin space and Qn{c) is the isospin coupling factor for the nucleon pole as shown 
in Table 2. 

(ii) Delta isobeir intermediate state and pion exchange 

^I(C) = -Qlic) (^) (^] ^(P3)757A?V(P1)^.(?) 

xi^{p^)p^.G%{pi)q^ij{p2), (21) 
9 



where G^j^{pi) is the delta isobar propagator as discussed in the previous subsec- 
tion, and Q\{c) is the isospin couphng factor for the delta pole as shown in Table 
2. 

(iii) Nucleon intermediate state and p meson exchange 



xV^(p4)757AP^^iv(A)(7a - :r^(^apq^MP2), (22) 

ZniN 



where Dp{q) is the propagator for the p meson defined as 

•^Mr - [-^] . (23) 

with g^^ being the usual metric tensor. In Eq. (22) a^^^ is defined as 

cr„v = -jinnlv-lul„), (24) 

Other quantities remain the same as earlier. 

(iv) Delta isobair intermediate state and p meson exchange 

xD'^/{q)^{p,)plG^MiQ''9''' - q'gnwi^iP^), (25) 

In Eq. (25) all the quantities are the same as defined earlier. 

(v) Nucleon intermediate state and a meson exchange 

AUc) = -QWa(—]i^iP3mpi)D.{qmP4) 

X^57t.P^DN{pi)llj{p2), (26) 

where D„{q) is the propagator for the sigma meson whose form is the same as 
that given by Eq. (19) except that in the denominator the pion mass is replaced 
by that of the sigma meson. 

(vi) Nucleon intermediate state and u meson exchange 
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In this case the form of the amphtude is the same as that of the p meson 
exchange. 

(vii) Pion emission from the decay of p meson in the intermediate 
states (Fig.lf) 

f - ik 

xDM^iP4H^r,q'MP2), (27) 

where q' = P4 — P2- 

(viii) Nucleon intermediate state and heavy ctxial meson exchange 

X ^ (Pi) Ib^fiPn^N {Pi)lbluiJ (P2) , (28) 
where is the propagator for the axial vector meson which is defined by 

D7{q) = -^f^^l (29) 

In Eq. (29), the mass of the axial meson was taken to be very large (188 
GeV) , as the corresponding amplitude is that of the contact term. 

The amplitudes given above can be simplified by contracting out the gamma 
matrices using the Dirac equation whenever applicable. The simplified expres- 
sions are given in Appendix A. The total amplitude is obtained by summing 
(coherently) the amphtudes corresponding to all the graphs. It should be noted 
that the exchange graphs have an extra minus sign due the antisymmetrisation. 

The general formula for the invariant cross sections of the N + N = N + 
N + n process is written as 
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where Afi represents the sum of all the amplitudes, Pj and Pf the sum of all 
the momenta in the intial and final states respectively, and pa the momenta of 
the three particles in the final state. The corresponding cross sections in the 
laboratory or center mass systems can be written from Eq. (30) by imposing the 
relevant conditions. 

3 Comparison to data and discussions 

The model presented in the previous section has been used to study the available 
data on the total, double and triple differential cross sections for the p + p ^ 
n + 7r+ + p and p+p— >p + 7r°+p reactions for beam energies 300 MeV to 2 
GeV. We emphasis here that the parameters described in the previous section 
have been kept fixed throughout in all the calculations described subsequently. 

3.1 Total pion production cross section 

In Figs. 2a and 2b, we show a comparison of our calculations with the exper- 
imental total cross sections for the reactions p + p— ^p + n + Tr"*" and 
p + p— >p-|-p-|-7r° respectively, as a function of beam energy. The dashed 
lines represent the results obtained by including only those graphs where pion 
production proceeds via intermediate delta isobar states, while the dashed-dotted 
lines give the results with only nucleon intermediate states. The solid lines shows 
the results where all the graphs are included. We note that the measured cross 
sections are reproduced reasonably well by our calculations in the entire range of 
beam energies. This is remarkable in view of the fact that none of the parameters 
of the model have been adjusted to the pion production data. Furthermore, they 
were determined by fitting the NN elastic scattering data (differential cross sec- 
tions) at beam energies above 1 GeV. The same set of parameters (with an energy 
dependence of the coupling constants given by Eq. (11) seems to reproduce the 
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data at low energies as well. 

The contribution of the delta isobar excitation is not important at beam ener- 
gies below approximately 350 MeV, which is due to the fact that at lower energies 
pions are predominantly in a relative S'-state; thus the possibility of forming a 
delta isobar is greatly reduced. However, the pion production is dominated by 
the A isobar excitation at higher beam energies. 

In Fig. 3, we show the contribution of various meson exchange processes to 
the total cross section of the {p,mT^) reaction in the considered range of beam 
energies. The contributions of the heavy axial meson exchange are not shown in 
this figure as they are negligibly small. We note that the pion exchange graphs 
dominate the production process for all the energies. The contribution of the 
p meson exchange is almost negligible at lower beam energies, and even in the 
energy range of 1 - 2 GeV, it is at least an order of magnitude smaller than 



that of the pion exchange process. Dmitriev et al. have made the similar 
observation for the reaction p + p n + within a similar type of model. 



This lends a posteriori credence to the calculations presented in Ref. |T0[ for the 
A{p, 7t)B reaction within a covariant two-nucleon model at a beam energy of 800 
MeV, where only one-pion exchange intermediate processes were considered. 

At the incident energies very close to the pion production threshold the con- 
tributions of the a and uj meson exchange are as strong as that of the pion 
exchange. Therefore calculations, which do not include these mesons, underesti- 
mate the pion production cross sections near the threshold (see eg. Horowitz et 
al. in [13], and 0). 

As we note from Table 1, the cut-off parameters A must have values in the 
range of 1.0 - 1.6 GeV in order to fit the NN scattering data. Similar values for this 



parameter have been used in almost all the earlier calculations ||T0|, |3^ , |33| , |3J] of 
the proton induced pion production (which mostly included only pion exchange 
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contributions, requiring therefore only A^r and Aa which were taken to be the 
same). In Fig. 4 we show the total cross sections for the p{p, mT~^)p reaction as a 
function of beam energy for two values of A^^. The solid and the long dashed lines 
show the results obtained with values 1.005 GeV and 0.63 GeV respectively for 
this parameter. It is apparent that the latter value of A^^ leads to a smaller cross 
section and a poorer fit to the data. Smaller values of A imply a reduction in the 
range of the pion exchange contribution to the NN interactions which reduces its 
contribution to the pion production cross section as this is the dominant process 
as has been discussed earlier. Although a value of A^ < 0.8 is consistent with 
the chiral bag model |3^, the quantitative description of the NN data requires a 



value > 1.0 GeV We therefore, stick to the value of A^^ as shown in Table 1. 

In Fig. 5, we investigate the effect of using the pseudoscalar coupling for 
the nNN vertex on the 7r+ production cross section. In this figure we show the 
ratio of the total cross sections obtained by using the pseudovector (apy) and 
pseudoscalar {crps) couplings for the ttNN vertex as a function of beam energy. 
The solid line represents the results obtained by including all the graphs while 
the long dashed line give the results obtained with only nucleon intermediate 
excitations. It is clear that the pion production is enhanced if the pseudoscalar 
coupling is used. This effect is very large at lower beam energies. At higher beam 
energies, where A isobar excitation dominates, this effect still persists even when 
all the graphs are included in the calculations. 

In the one pion exchange picture of the NN interaction, when both nucleon 
are on the mass shell, it will not be possible to make a difference between the 
pseudovector and pseudoscalar couplings, as both will produce the same result. 
However, in Fig. 1, one of the nucleons can go off-shell. Therefore, in this case 
the two couplings will produce different results. One can show (taking the case 
of just one diagram, say for example Fig. la) that the ratio of the amplitudes 
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obtained by using pseudovector and psedoscalar couplings is given approximately 
by (i?^/2m7v). At lower beam energies this ratio is very small, which explains to 
some the extent the trend seen in Fig. 5. 

In general, the pseudovector coupling for the irNN vertex is preferred. The 
p-wave interaction between pions and nucleons is very strong and attractive; 
the formation of the A isobar is the consequence of this interaction. On the 
other hand, the s-wave interaction is very weak. The pseudovector coupling 
automatically incorporates these features due to the derivative term. One can 
show purely on formal grounds that it is also consistent with the constraint 
imposed by the PCAC hypothesis [|l^]. The derivative of the axial vector current 
operator based on the pseudovector ttNN coupling vanishes as the pion mass 
goes to zero. The pseudoscalar coupling does not satisfy this criterion; it also 
produces a s-wave interaction which is much too strong [0. Furthermore, this 
coupling has a history of problems in Dirac calculations [^]. Therefore we do 



not consider the pseudoscalar couplings for the ttNN vertex in our calculations. 
3.2 Double differential cross section 

We have also calculated the double differential cross sections for the pp tt+X 
reactions and have compared the results with one set of existing experimental 
data at the proton incident energy of 800 MeV, in Fig. 6; this set of data 
were also analysed by Verwest |T5[ . It is clear that the contributions of the delta 



isobar excitation (dashed lines) dominate the total cross sections (solid lines), 
whereas those of the nucleon intermediate state (dotted lines) are very small. 
The agreement between theory and the data is reasonably good at forward pion 
angles. However, at higher pion angles our calculations overpredict the data 
by factors of 1.5 to 2. The peak cross sections are affected by the different 
multipolarities (dipole or monopole) of the form-factors (Eq. (10) - (11)) and 
forms of the delta decay width (Eq. (15)). For instance, with a monopole 
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form factor for the the NAp vertex, the cross sections near the peak are reduced 
by approximately 20 %. On the other hand, using the width r^^er enhances the 
cross sections near the peak region further by about 5-6 %. There is, therefore, 
some scope to explain the small disagreement between data and out calculations 
at larger pion angles. This also indicates that these measurments may be useful in 
differentiating between various forms of the formfactors and delta decay widths. 

The experimental points towards the larger momentum ends of the spectra 
come from the processes like nd final states which is obviously not included in 
our calculations. 

3.3 Triple differential cross section 

We now compare our results with data taken in kinematically complete measur- 
ments where all the three particles in the final channels are measured. In Fig. 
7, we present a comparison of our calculations with the data p for the triple 
differential cross sections for two sets of pion and proton angles as a function of 
outgoing proton momentum at the beam energy of 800 MeV. The peak reagion 
in this figure corresponds to the final state proton momenta associated with the 
formation of the delta isobar and its decay. The solid lines (dashed lines) in 
this figure show the results of calculations which include the contributions of all 
the graphs (only delta isobar intermediate states). Clearly the peak region is 
dominated by the delta excitation process. At higher angles and in the regions 
away from the delta excitation peak, the contributions of nucleon excitation are 
not negligible as can be seen in the lower part of this figure. The shift in the 
peak position towards higher momenta at larger angles is due to the fact that at 
larger pion angles the delta isobar formation is associated with larger outgoing 
pion momenta. 

It should, however, be noted that the position of the peak in the calculated 
cross sections is shifted towards somewhat larger momenta as compared to the 
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experimental data. Obviously the form of the delta decay width will play a crucial 
role in determining the peak position. In Fig. 7 we have used the decay width 
r^er- The effect of using other forms of the deacay width is shown in Fig. 8, where 
the solid lines represent the results of the calculations performed with the delta 
decay width of Verwest, while the dashed lines corresponds to that of Dmitriev. 
The short dashed lines show the results obtained with Dmitriev's form but with a 
value of Fo = 0.100 GeV [|]. We see that the peak positions of the experimental 
cross sections are correctly reproduced by calculations using the Decay width of 
Dmitriev. Hoewever, the absolute magnitudes of the cross sections are a bit too 
high at larger pion angles. 

4 Summary and conclusions 

In this paper we have presented a fully covariant one-boson exchange model 
to describe the p{p,nn^)p and p{p,pTT^)p reactions in the beam energy {Eheam) 
range of 300 MeV to 2 GeV. The model contains only the physical parameters (like 
coupling constants and cutoff masses), which were determined by fiting to the NN 
elastic scattering data at three beam energies above 1 GeV. All the parameters so 
determined were held fixed throughout the considered energy range; none of the 
parameters were adjusted to the pion production data of any kind. We considered 
the exchange of the vr, p, a and oj mesons; the latter two were not considered 
in earlier calculations of the NNtt processes. The calculations included pre- and 
post-emission graphs and considered the excitation of both the delta isobar and 
nucleon intermediate states. The deacy of the p meson in flight (the so-called 
intermediate emission diagram) was also incorporated in our calculations. 

We found that the pion exchange processes dominate the cross sections in 
the entire energy region. The contribution of the p meson exchange is almost 
negligible at lower beam energies. Even in the energy range of 1 - 2 GeV its 
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contribution is less than 10 %. This is an important observation as it imphes 
that in the description of the A{p, tt)B reaction at beam energies around 1 GeV 
within the covariant two nucleon model, it suffices to consider only pion exchange 
intermediate states, which reduces the intricacy of the calculations quite a bit 
^. The exchange of a and u mesons is important for beam energies closer to 



the pion production threshold, while at intermediate and higher beam energies 
their contributions are of the same order of magnitude as that of the p meson 
exchange. Therefore, it would be important to include the exchange of both cr 
and uj mesons in order to explain the near threshold data on pion production in 
pp collisions taken recently at the Bloomington cooler cyclotron f^. It should 
be noted that a fully covariant calculation is necessary at even these energies 
because the non-relativistic reduction of the ttNN Lagrangian has ambiguities, 
which cast doubt on the results obtained within such approaches [p!0|, 

The excitation and decay of the delta isobar dominates the pion production 
processes at Ebeam > 0.5 GeV where it accounts for almost entire cross sections 
on its own. However, this process is not important for lower beam energies; it 
contributes almost negligibly to pion production near the threshold. 

We have also reviewed the derivation of the covariant delta isobar propagator 
and have showed that the propagator given by Williams does not represent the 
correct propagator for a massive spin-3/2 field. This point was further stressed 
by performing numerical calculations with this propagator. We have found that 
the Williams propagator leads to very large cross sections ( which look almost 
like a pole ) for total as well as differential cross sections in certain energy regions 
which is clearly unphysical. 

The calculations performed with a pseudoscalar coupling for the nNN ver- 
tex were found to produce too large cross sections particularly near the pion 
production threshold. This is due to the fact that pseudoscalar coupling alone 
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(without any crvr coupling) produces a much too large S'-wave nN interaction, 
which dominates the pion production near the threshold. However, this is clearly 
not in agreement with the experimental data. The pseudoscalar coupling is not 
consistent with the constraints of the PCAC. Therefore, we have preferred to use 
the pseudovector couplings in our calculations. 

We have also analysed the available double and triple differential cross sections 
for pion production and found that our model provides a reasonable description 
of these data as well. The triple differential cross sections are sensitive to the 
form of the momentum dependent delta decay widths and probably also to the 
form of the cut-off parameters. As there are many versions of this width available 
in the literature, the complete kinematical detection of three particles in the final 
channel may perhaps help in making distinction between them. This would lead 



to a better understanding of the off-shell behaviour of the delta isobar [pj] . 

In this work we established a covariant framework to describe the inelastic 
channels of the NN collisions, and used this to investigated the strongest of them, 
the pion emission. This should now form the basis for the study of relatively 
weaker inelastic processes like emission of t] mesons and also the rj' meson which 
is thought to provide a novel tool to study the baryonic resonances around 2 GeV 

m 
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Appendix A 

By using the complicated but straight-forward algebra of 7 matrices and the 
Dirac equation for a free particle the amplitudes given in subsection (2.3) can 
be rewritten into relatively simple forms which are suitable for numerical calcu- 
lations. We have dropped the isospin factors in the expressions shown below, 
however in actual calculations they are included. 

(i) Nucleon intermediate state and the pion exchange 

AUc) = -( — ] 2 ^ 2 2 \ ^(P3)75^(Pl) 

x^{p,){E + F^^p';^)^{p2), (31) 

where 

E = 2mjv(m^-p2) (32) 

F = p'i+3m%. (33) 

It may be noted that if the pseudoscalar coupling for the ttA^A^ vertex is used, 
we get the same expression for this amplitude but the constants E and F are 
defined in the following way 

E ^ (34) 
F = 4:171%. (35) 

(ii) Delta isobeir intermediate state and the pion exchange 

f\ff*V 1 1 



xHPsHi^{PiMP4){G + i/7X)V'(P2). (36) 

In Eq. (36) the expressions of G and H will depend on the form of the delta 
isobar propagator used in the derivations. For the propagator (Eq. (14)), 
one gets 
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G = ruNqp^ - -niNPiP^ + -niNm^ — ^ — vhn + — 

PiPnml piqml 2 1 ^ 
^ \- ^ \- ruAqPn - i^mAPiPn + i^ruAm^ 

"^PiP^PiQ 
3mA ' 

„ 2 2 ,1 2 '^PiPirPiQ PPniriN 

PiQUlN 2 1 1 

— o ^ ornATTlN + IzPiPn - l^PiQ- 

6mA <J 6 6 

Whereas with the propagator (see appendix B), we obtain 

ri 2 1 2 '^PiPi^PiO. , 2(pp^)2 

G = vfiNqp-K - -mNPiPTT + -^mNm^ — niN H Tr^ruA 

3 3 Spf Spf 

PiP^ml Piqml 2 1 ^ 

rriA + „ 2 + mAqPw - ^mAPiPT, + ^mAm^ 

6pi 3pf 3 3 

"ipiPnPiq 

„ 2 2 1 2 '^PiP-KPiQ PpT^rriN 

H = gp. - -m, + -m, - ^-^ ^mA 

PiqmN 2 1 1 

^-^"^A + ^mATriN + -PiP^ - -Piq. 

6pf 3 3 3 

(iii) Nucleon excitation and p meson exchange 
We make use of the following relations 

4*- -9^1 - ». 



^{P3h^^{Pl){P3 - Pi) = 0, 



to reduce the amplitude given by Eq. (22) as follows 

Kic) = 9l f — ) 2 ^ 2 2 \ V^(P4)(i^l - ^2)V'(P1). 

" \m^J r - mf^pf - my 
In Eq. (44) Di and D2 are defined as 



^1 = i^fiPn^uPi^v^'' - "^iV7/.P^^7r,^'')75, 
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where 



-2pip^-fi^b^-f„q''j5 + mN^i^p'i7ub''-fr,q''75] , (46) 



k 

b^, = -0(^3) [(1 - kp)-f^ + — (p3m - qpmMiP^)- (47) 



(iv) Delta isobar excitation and p meson exchange 

Using Eqs. (41) - (43) we can write the amphtude given by Eq. (25) as 

A%{C) = -,,^5^(^)^(p3)[(l-Mr + -^P3l^(pO 

^^{pAW.G%^,q'l, - PlG^q'l.l,). (48) 

This equation can be further simphfied by using the forms of the propagator G% 
as discussed earlier. These expressions are not being given here as they are very 
lengthy even though it is straightforward to derive them. 

(v) Nucleon excitation and a meson exchange 

\m^J q^-mipt-m% 

X'0(p4)(2mJv7MP$^ + 1^- 2pjP7r)75V'(P2)- (49) 

(vi) Pion emission from the decay of p meson in the intermediate 
states 

AP^"" = i2mNgpgp^^ I —] — — o 

\mT, J — nip g'^ — 

x'0(P3)[(i - ^p)7MP^^ + —{psPn — r?^^T)]'0(pi) 

niN q 

X'0(P4)75^(P2)- (50) 

(vii) Nucleon intermediate state and the exchange of heavy ctxial 
meson 
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y ) — m\ pf — ml 
xi^ipdYb^ - 2m^tlj{pi)YPnuYbMP2), (51) 

where 

= '^(Pshblt^^iPi) (52) 

The Numerical evaluation of the amplitudes written in the forms given above can 
be carried out very effectively by using the techniques of the Clifford algebra. 

Appendix B 

In this appendix we present the discussion on the form of the delta isobar prop- 
agator. 



The free Lagrangian density for the massive spin-3/2 field is written as ||3] 

Ca = ^'^A;..^^ (53) 
where the most general form of A^^, is given by 

= -[{-idf.-f^' + mA)gf,u-'i'Zi{-ff,d^ + df,-f^) (54) 
-Z27^((9a7^)7,. - Z^mA'jf.'ju] 

Here Zi is an arbitrary parameter subject to the restriction that Zi ^ —1/2, 
and Z2 and Z3 are defined as 

Z2 = ^{SZf + 2Zi + 1), Z3 = (SZl + 3Zi + 1) (55) 

Physical properties of the free field do not depend on the parameter Zi, which 
is chosen to be real. This is due to the fact that £a is invariant under the point 
transformation H3I 



^ ^'^ + 07^7'^^^ (56) 
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where a 7^ — |, but is otherwise arbitrary. 

The wave equation for the spin-3/2 particle (the Rarita-Schwinger) equation 
is written as 

A^,(p)^'^ = (58) 

Operating on Eq. (58) with 7^ and d'^ we get the local wave equation for a 
spin-3/2 particle along with the constraint equations 

{id^-rriAW = (59) 
7/.^'' = (60) 
d^"^^ = (61) 

It should be noted that Eqs. (59) - (61) are obtained only when the restriction 
Zi 7^ —1/2 and the definitions of Z2 and Z3 as given in Eq. (55) are used. 

The propagator for a massive spin-3/2 particle satisfies the following equation 
in the momentum space 

^f^Ap)GAa = -9f^a (62) 

Solving for G and making the particular choice of Zi = — 1 we get the following 
equation for the delta isobar propagator 

G^ip) = _ ^(^ + "^^) [^ _ 1^^^ _ -^p + -l^(p 7^ - p^7 )] (63) 

A slightly different form of the A propagatator has been suggested by Williams 
2^, which is given by 



„AM// X i(-/> + mA), 12 1 „, ,T . ^, 

= ~ p2_^2^ t^M'- - 3^m7z. - ^PmP- + ^Pr,7^(PM7<. - P.7/.)] (64) 



Williams has argued that the Eq. (64) should be the valid form of the A isobar 
propagator both on and off the mass shell. However, we show in the following 
that it can not be the correct spin-3/2 propagator. 
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First, we note that if is the propagator for a spin-3/2 particle, then it 
should satisfy 

M^'''{p)GZ^^ = (65) 

where A' is defined as 

A^^ = -[{-id^Y + ^A)g^^,y + ^>^l^,^„-iX^^-f^], (66) 

with the limit that A ^ oo (known as Feynman gauge). However, Eq. (66) is not 
consistent with the general form of the Lagrange function of a spin-3/2 particle 
as given by Eqs. (54) and (55). 

Second, in Ref. [^, it has been shown that the propagator has no 



inverse, and thus it can not be the propagator of a spin-3/2 particle. A similar 



difficulty is associated with the propagator suggested by Adelseck et al. ||21| . 
In Fig. 9a and 9b, we study the effect of using the two delta isobar propagators 
and G^^ in the calculations of the total and triple differential cross sections 
for the p{p, n7T~^)p reaction to stress further the difficulties associated with the 
Williams propagator. The solid lines show the results obtained by using the 
former propagator while the dashed line the later one. It is clear that the G^J^ 
produces very large cross sections in a certain range of beam energies due to 
the fact that the factor p"^ present in the demoninators of some terms of this 
propagator, becomes very small. This effect is particularly very strong for the 
pre-emission graphs. For example, for such a graph which is analogous to Fig. 
la we have 

2 / 2 2\ 

Pi = {P2-PJ 

= m% - 2E^mN + ml (67) 

It can be seen that for E^^ ^ 0.479 GeV , Pi is very close to zero. This energy 
corresponds to an incident proton energy of approximatly of 1.2 GeV. Therefore, 
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at incident energies around this value the terms which are proportional to 
in become very large in the pre-emision graphs which in turn leads to 

huge total cross sections. In the post-emission graphs this problem is not so 
severe. This very clearly shows the difficulty that one encounters while using the 
propagator . It must be mentioned here that normally the contributions of 
pre-emission graphs are much smaller than those of the post-emission ones if one 
uses the propagator G^^. 

The effect of using G^^ is very drastic in case of triple differential cross 
sections as can be seen in Fig. 9b. Williams propagator leads to a very large 
cross section (which looks almost hke a pole) at certain value of the outgoing 
pion momenta due to the same reason as discussed above. 
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Tabic 1: Coupling constants for the NN meson vertices used in the cal- 
culations 



Meson 




e 


A 


mass 








(GcV ) 


(GoV) 


TT 


12.562 


0.1133 


1.005 


0.138 


(7 


2.340 


0.1070 


1.952 


0.550 


UJ 


46.035 


0.0985 


0.984 


0.783 


P 


0.317 


0.1800 


1.607 


0.770 


kp = 6.033, = 0.0, g,^^ = 2gp 
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Table 2: Isospinfactors for pole diagramms 



nucleon pole 




delta pole 


pp — >• np7T~^ 


PP — >• npTT"*" 


graph 


isovector 


isoscalar 


graph 


isovector 


isoscalar 


a 


a/2 


a/2 




a 


-a/2/3 





b 










b 


V2 





c 










c 


V2 





d 


V2 


V2 




d 


-V2/3 





pre a 


-V2 


V2 




pre a 







pre b 


2^2 







pre b 


72/3 





pre c 


2^ 







pre c 


V2/3 





pre d 




V2 




pre d 


V2/3 





PP ppTl° 




PP —>■ ppn° 


all graphs 


1 


all graphs 


2/3 



Table 3: Isospinfactor for direct diagramms 



intermediate 


PP 


npn^ 


graph 




e 


-iV2 


f 


iV2 


g 


iV2 


h 




PP — 


ppiT° 


all graphs 
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Figure Captions 

Fig. 1 . Feynman diagrams for emission of a pion in the nucleon - nucleon col- 
lisions, (a) Exchanged meson starts from nucleon 2 (momentum P2), is 
absorbed by nucleon 1 (momentum pi) which is excited to the nucleon or 
delta isobar intermediate state with momentum {pi) which then decays into 
the nucleon 3 (momentum p^) and the outgoing pion (momentum p^). The 
nucleon 2 goes on to nucleon 4 with momentum (^4). (b) Exchange part 
of the diagram (a), (c) Same as (a) but the exchanged meson starts from 
nucleon 1 and is captured by nucleon 2 which is excited to intermediate 
states which then decay into nucleon 4 and the outgoing pion. Nucleon 1 
goes on to become nucleon 3. (d) Exchange part of diagramme (c). (e)-(h) 
Direct and exchange diagrams showing the processes where the exchanged 
p meson starting from one of the interacting nucleons decays into two pions 
in flight. One of them is the outgoing pion, the another one is absorbed 
by the other nucleon. These diagrams are referred as intermediate graphs 
in the text. Note that there are also the pre-emission counter-parts of the 
diagramms (a) - (d) where pion are emitted before collisions. These graphs 
are not shown here, but their contributions are included in the calculations. 

Fig. 2a . The total cross section for the p{p, n'K^)p reactions as a function of beam 
energy. The dashed and dashed-dotted lines represent the results of the 
calculations performed with only delta isobar intermediate states and only 
nucleon intermediate states, respectively. The sum of all the graphs is 
represented by the solid curve. The experimental data are taken from 0. 

Fig2b . Same as Fig. 2a but for the reaction p{p,pTT^)p. 

Fig. 3 . Contributions of various exchanged mesons to the total cross section for 
the reaction p{p,p'iT~^)p as a function of beam energy. The long dashed, 
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dotted, short dashed and dashed-dotted curves represent the contributions 
of TT exchange, p exchange, uj exchange and cr exchange, respectively. 

Fig. 4 . Total cross sections for the reaction p{p,pn'^)p as a function of beam 
energy for two values of the cut-off parameter of the ttNN vertex. The 
solid (long dashed) line is the result of calculations performed with a value 
of 1.005 GeV (0.631 GeV) for this parameter. 

Fig. 5 . Ratio of the total cross sections calculated with pseudovector and pseu- 
doscalar couplings for the ttNN vertex for the same reaction as in Fig. 4, 
as a function of beam energy. The long dashed line represents the results 
obtained with only nucleon intermediate states while the solid line contain 
all the graphs. 

Fig. 6 . The double differential cross section for the reaction pp ^ n'^ + X as 
a function of pion momentum for pion angles of 20°, 40" and 60" at a 
beam energy of 800 MeV. The dashed (dotted) hne represents the results 
obtained with only delta isobar (only nucleon) intermediate states. The 
solid line shows the results obtained by including all the graphs. 

Fig. 7 . The triple differential cross sections for the p(p, nvr^) reaction at the 
beam energy of 800 MeV as a function of the outgoing pion momentum. 
The upper part shows the results for the proton and pion angles of 15° 
and 21° respectively while the lower part for 25° and 40° respectively. The 
dashed lines show the results when only delta isobar intermediate states are 
included in the calculations while the solid lines depict the results obtained 
by including all the graphs. 

Fig. 8 . The effect of using various forms of delta decay width in the calculation 
of triple differential cross sections for the same reaction as in Fig. 7 and 
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at the same beam energy. The sohd and dashed hnes represent the results 
obtained with the forms of the delta decay widths given by Verwest and 
Dmitriev with a free delta decay widths of 120 MeV. The dotted lines 
represent the results obtained with the later but with a a free delta decay 
width of 100 MeV. The results are shown as a function of outgoing pion 
momentum. 

Fig. 9a . Total cross section for the p(p,p7r+)p reaction as a function of beam energy 
calculated with the propagators of Benmerrouche et al. ||2^ (solid line) and 
Wilhams pO| (long dashed line). 

Fig. 9b . Triple differential cross section for the p{p,pTT~^)p reaction as a function 
of outgoing proton momentum corresponding to proton and pion angles of 
15° and 21° respectively, The solid and long dashed curves have the same 
meaning as in Fig. 9a. 
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